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Abstract 



The BATSE experiment on GRO has demonstrated the isotropic arrival directions 
and flat log N vs. log S of cosmic gamma-ray bursts. These data are best explained if the 
burst sources are distributed throughout an extended spherical Galactic halo, as previously 
suggested by Jennings. The halo's radius is at least 40 Kpc, and probably is more than 100 
Kpc. I consider possible origins of this halo, including primordial formation and neutron 
stars recoiling from their birthplaces in the Galactic disc. A simple geometrical model 
leads to a predicted relation between the dipole and quadrupole anisotropy. I suggest 
that neutron stars born with low recoil become millisecond pulsars, while those born with 
high recoil become the sources of gamma-ray bursts; these populations are nearly disjoint. 
Quiescent counterparts of gamma-ray bursts are predicted to be undetectably faint. 
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The first results from the BATSE on GRO (BATSE Science Team 1991) have revived 
the question of the distribution of gamma-ray burst sources in space. Their chief results, 
isotropy of gamma-ray burst directions and a log N vs. log S slope significantly flatter than 
-1.5, confirm earlier reports (see, for example, Meegan, Fishman and Wilson 1985 and the 
review by Cline 1984). Questions of relative calibration of different instruments and the 
paucity of good directional data permitted skepticism in the past. Such skepticism is no 
longer tenable, and the theoretical questions raised earlier must be faced. 

An isotropic distribution of sources implies that, out to the maximum distance of 
observation permitted by instrumental sensitivity, all directions contain equivalent source 
populations. The source population for an observed flux or fluence S is expressed as the 
integral 

N(Cl,S) = J n(n,r,4%r 2 S)4Tvr 2 dr, (1) 

where n(0, r, 47rr 2 S) is the density of sources at a distance r from the observer in the 
direction O radiating a total power or energy 4nr 2 S; a temporal average is implicit, and 
it is assumed that the sources radiate isotropically. While it is possible to obtain N 
independent of O for a variety of artificial spatial distributions n, the only plausible n 
which does not require the observer to be in a preferred position (for example, at the 
center of spherical shells) is a homogeneous distribution out to the greatest distances at 
which the most luminous burst may be observed. Such a distribution implies N oc S~ 3 / 2 , in 
contradiction to observation. This dilemma was fully appreciated and discussed by several 
of the contributors to the 1983 Symposium on High Energy Transients in Astrophysics 
(Woosley 1984). 

The only geometric resolution of this problem is the Ptolemaic gamma-ray burst 
universe, as discussed by Cline (1984). The flat dependence of N on S (equivalent to 
a peak in a V/V max distribution at 7 C V max , and to a deficiency of faint bursts) is 
explained if the sources are distributed throughout a volume of finite extent, with most 
bursts within the volume luminous enough to be observed. This volume is spherical, with 
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the Earth at its center. The sources must be at cosmological distances, or be distributed 
in an extended halo around our Galaxy (such as the Massive Extended Halo suggested 
by Jennings 1984), or be evidence of a previously unsuspected Earth-centered spherical 
distribution. 

At cosmological distances a typical burst of 10~ 7 erg/cm 2 fluence must have radiated 
~ 10 50 erg in gamma-rays. This is energetically consistent with the spiraling-in or coales- 
cence of binary neutron stars, if they convert gravitational energy to gamma-rays, perhaps 
through a magnetohydrodynamic dynamo, or with a relativistic supernova shock breakout, 
but the mechanisms are obscure. The breakout of a relativistic shock from a spherical shell 
of radius R s would produce a pulse of radiation whose envelope would be proportional to 
the local (relativistically collimated) emission pattern f(Q(t)), where 9{t) = (2ct/R s ) 1 ^ 2 
and #~7 _1 ^Clis assumed. The rapid fluctuations of intensity under a fiat envelope 
observed in many bursts might then be explained by shock refraction in a turbulent stellar 
envelope. Unfortunately, the implied values of R s ~ 2d^7 2 , where tb is the burst duration, 
are unreasonably large, particularly because the 7-7 pair-production constraint on 100 
MeV gamma-rays implies 6 < 0.01. 

Such catastrophic events might be expected all to have similar time histories, or to fall 
into a small number of homogeneous classes, inconsistent with observation. The observed 
soft gamma repeaters, the 8 second periodicity of the March 5, 1979 burst, the absence 
of bursts from local supernovae, the identification of 400-500 KeV spectral features with 
positron annihilation lines of modest redshift, and the complex and varied time histories 
of intensity within bursts all disfavor these cosmological hypotheses. 

If an observer is displaced by a distance a from the center of a sphere of radius 
R^> a, the mean cosine of the angle 9 between the points in the sphere and the direction 
of displacement is found by integration, keeping only the lowest order terms in a/R: 



where 6 = 6>o + A6> 6>o — 7 sin and the integration runs over coordinates (r, 6>o,</>) 
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measured from the center of the sphere. Adopting a = 8.5 Kpc for our displacement from 
the center of the Galaxy and (cos#) < 0.08 (BATSE 1991; lex bound) implies R > 106 
Kpc; the 3cr bound of (cos 6) < 0.22 relaxes this to R > 39 Kpc. There is little other 
direct evidence for a spherical component of our Galaxy extending to these radii, but it is 
believed (based on dynamical measurements) that the mass distributions of many galaxies 
similarly extend far beyond their visible radii. The requirement that the halo be spherical 
constrains models of its formation and gravitational potential. 

It is possible that the objects (by consensus, probably neutron stars) producing 
gamma-ray bursts were formed in this extended halo early in the history of the Galaxy, 
and that they remain active for > 10 10 years. The total mean gamma-ray burst flux cross- 
ing the Galactic plane is ~ 10 -11 erg/cm 2 sec, corresponding to a Galactic luminosity of 
~ 10 37 erg/sec even for a very extended halo (R ~ 100 Kpc), or to ~ 10 26 erg/sec M & for 
a halo mass of ~ 10 11 M Q , comparable to the mass of the visible Galaxy. The cumulative 
power radiated is ~ 3 x 10 43 erg/M©. If all the mass were in neutron stars, this energy 
would be comparable to their magnetospheric energy for B ~ 10 13 gauss, or to their ro- 
tational energy with moments of inertia of 10 45 gm cm 2 and spin periods ~ 25 seconds. 
Presumably only a fraction of the mass is in neutron stars; their expected magnetospheric 
energy is inadequate, but their rotational energy may be sufficient. 

Alternatively, the sources of gamma-ray bursts may be born within the visible Galaxy 
but escape to an extended halo. It is well known (Lyne, Anderson and Salter 1982) that 
many (but probably not all) pulsars are produced with recoil velocities ~ 200 km/sec, 
sufficient to escape the Galactic disc. When combined with their initial orbital kinetic 
energy, they may be sufficiently energetic to escape into an extended halo. The time 
required to fill a halo of radius R ~ 40-100 Kpc is 1-3 xlO 8 years. This hypothesis is 
complicated by the fact that their initial orbital velocities about the Galactic center will 
tend to concentrate the neutron stars to low Galactic latitudes, and any low-recoil fraction 
will remain in the Galactic plane, concentrated toward the Galactic center. 
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The lifetimes of these neutron stars as gamma-ray burst sources must not be less 
than this filling time (or the bursts would be concentrated within the visible Galaxy, and 
peak towards its center) , but could be arbitrarily greater. It is not possible to estimate the 
number of neutron stars required because their burst frequency is unknown. If the neutron 
stars freely escape to infinity and radiate indefinitely their n oc r -2 implies an asymptotic 
d log N/d log S = —0.5, not as flat as the dlogN/dlogS = found for a confined space 
distribution, but flat enough to explain the data. The energy required is the same as in 
the case of gamma-ray burst sources born in the halo. 

For a simple model of gamma-ray burst sources recoiling from their birthplaces in the 
Galactic disc, it is possible to calculate the expected dipole moment of their distribution 
of arrival directions with respect to the Galactic center and the quadrupole moment with 
respect to the Galactic plane. Assume that we reside a distance ro from the center of a 
Galactic disc of radius L, uniformly populated with source birthplaces, and that each point 
in the disc uniformly fills with observable gamma-ray bursts a sphere of radius R ^> L, ro 
with itself at its center. The uniformly filled sphere is a rough approximation to the motion 
of neutron stars in the unknown potential of the Galactic halo. Using the result (2), which 
now describes the observed asymmetry of each sphere, with a 2 = r 2 + r 2 — 2rr cos <j> 
and (r, <f>) denoting the coordinates of the disc point, projecting this asymmetry onto the 
direction to the Galactic center, and integrating over the disc yields 

(cob<9)«^, (3) 

where 9 is the angle between a burst's arrival direction and that to the Galactic center. 
Application to the BATSE data yields the same bounds on R as discussed above, although 
now R refers to the radius of a sphere populated by gamma-ray bursts from a given 
neutron star birthplace, rather than the radius of a sphere of source positions centered 
at the Galactic center; the distinction is slight for R ^> L, ro, as must be the case. If 
the assumption of uniform filling is replaced by d~ 2 density distribution, where d is the 
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distance from the birthplace (reflecting an outward streaming distribution, rather than a 
trapped one), then the result (3) is multiplied by a factor O(ln(R/r )) ~ 2. 

It is similarly possible to calculate an average (sin 2 /?), where Q is a burst's Galactic 
latitude. For an observer in the Galactic equatorial plane but at a distance a from the 
center of a uniformly populated sphere of radius R, 

(sm z 3) « -, (4) 

x M/ 3 5R 2 K 1 

to lowest non- vanishing order in a/R. Averaging over spheres centered in the Galactic disc 
yields 

/•2m 1 4r 2 + 2L 2 

(sin ^ri^- (5) 

Because of the quadratic dependence on R the resulting bounds are less strict than those 
obtained from (2) or (3). 

The results (3) and (5) may be combined to predict a relation between these two 
angular averages, if the relation between L and r$ is known. Take, for simplicity, the 
reasonable estimate L = y/2r . Then 

(sin 2 /3)^-j|(cos#> 2 . (6) 

This relation offers a test of the simple geometrical model used here. This test is readily 
satisfied by the extant BATSE data, but may become more discriminating as more data 
accumulate. 

The gamma-ray burst of March 5, 1979 may, in part, be explained by the extended 
halo hypothesis. If most other bursts are at distances less than the thickness of the Galactic 
disc, as was often supposed before the BATSE data became available, then the March 5, 
1979 event was extraordinary, because its luminosity must have been ~ 10 5 — 10 6 times 
greater than that of a typical burst. It is hard to reconcile this interpretation with the 
absence of bursts similar to that of March 5, 1979, but occurring in our own Galaxy, which 
contains ~ 10 times the mass of the LMC, and might therefore be expected to have ~ 10 
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times as many such events (but closer and even brighter!). If, instead, bursts in general 
are distributed throughout a very extended halo of our Galaxy, extending to distances at 
least as great as that of the LMC, then the March 5, 1979 burst would be explicable as 
an unusually (but not extraordinarily) luminous member of a class of events with a broad 
distribution of luminosities. 

Another consequence of the hypothesis of the birth of gamma-ray burst sources in 
the Galactic disc is a further development of the empirical division of neutron star births 
into high and low recoil classes (Katz 1975). Millisecond pulsars, found in superabundance 
in dynamically fragile globular clusters (like their probable ancestors, the X-ray binaries) 
must have been formed recoillessly. Slowly spinning pulsars, most of which are single and 
are observed (Lyne, Anderson and Salter 1982) to have high space velocities (and hence 
high recoil) are the ancestors of gamma-ray burst sources. 

These categories are nearly disjoint because low-recoil millisecond pulsars retain their 
companions, initially suppressing their pulsar emission (Shvartsman 1971), but later spin- 
ning them up. The observed small magnetic dipole moment fx of millisecond pulsars reduces 
the power < (U 2 c/(3i?^ s ) potentially released by magnetospheric dissipation in a gamma- 
ray burst, explaining the absence of this Galactocentrically-concentrated population from 
the spatial distribution of bursts, although their large rotational energy suggests that a 
significant luminosity in microbursts may be possible. 

High-recoil neutron stars, born with slow spin but which lose any binary companion, 
must retain their magnetic fields for a very long time in order to produce gamma-ray bursts 
(empirically, Her X-l has retained its large field for at least 3 x 10 7 years, the time required 
for it to reach its height above the Galactic plane, and perhaps much longer). Low-recoil 
single neutron stars and high-recoil neutron stars which retain their companions, although 
individually interesting (Her X-l, PSR1913+16), are uncommon. 

If gamma-ray bursts are distant then they must be very powerful. At 100 Kpc a strong 
10~ 5 erg/cm 2 burst, radiating isotropically, emits 10 43 erg, posing energetic problems 
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similar to those posed by the March 5, 1979 event (Katz 1982). Nonthermal mechanisms 
involving the flare-like dissipation of much of a magnetosphere's energy are an attractive 
explanation of much gamma-ray burst phenomenology, but they are barely energetically 
adequate for a single burst. The possibility of regeneration of magnetospheric energy even 
during a burst (perhaps by twisting up the magnetosphere by differential rotation) must 
be considered. Note, though, that even if the burst population fills a sphere of radius R 
the brighter bursts will, on average, be significantly closer than R, reducing their energy 
requirements but increasing the spatial anisotropy of their distribution. 

A final consequence of the distances suggested here is that quiescent burst counterparts 
will be too faint to observe. Stellar companions (unless close and compact) will have been 
lost when recoil was imparted. Accurate burst coordinates and deep searches will show 
only blank fields, as has been the case for the few accurate positions presently available 
(Schaefer 1990). An essentially distance-independent argument (Katz 1985) implies that 
Schaefer's (1990; but see Zytkow 1990) reported optical flashes can only be explained as 
nonthermal emission. At a distance of 100 Kpc thermal reprocessing would be inadequate 
to explain the reported intensities (the inferred absolute magnitudes would be brighter 
than -15!). Similarly, at these distances self-absorption suppresses the reprocessed X-ray 
flux, explaining the low X-ray luminosity of gamma-ray bursts (Imamura and Epstein 
1987). 

I thank G. J. Fishman for supplying data in advance of publication and F. J. Dyson 
for discussions. 
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